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0^ We solve a nonlocal generalisation of the NJL model in the Hartree ap- 

proximation. This model has a separable interaction, as suggested by instan- 

ton models of the QCD vacuum. The choice of form factor in this interaction 

O , 

^Sj , is motivated by the confining nature of the vacuum. A conserved axial cur- 



X 



rent is constructed in the chiral limit of the model and the pion properties 
are shown to satisfy the Gell-Mann-Oakes-Renner relation. For reasonable 



> 

cn 
en , 

t^^ ' values of the parameters the model exhibits quark confinement. 

o' 

^ I. INTRODUCTION 



The Nambu-Jona-Lasinio (NJL) model |]l| is often used as a model to study the hidden 
chiral symmetry of QCD [Q. It describes fermions interacting via a local, chirally invariant 
four-point coupling. For large enough coupling strengths it leads to a vacuum state with a 
quark condensate and pions which are approximate Goldstone bosons. The local nature of 
the interaction considerably simplifies the calculations compared with ones based on finite- 
range forces. However it also means that loop diagrams are divergent and must be regularised 
in some way. 

This regularisation scheme forms part of the model and introduces a cut-off as an extra 
parameter. Many procedures have been used in the literature, for example: three- and 
four- momentum cut-offs, Pauli-Villars, proper-time regularisation [Q. All of these yield 
qualitatively similar results. None however has any clear interpretation in terms of the 
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underlying theory, QCD. Cutting off the quark propagator at high momenta, as most of 
these schemes do, does violence to the Hilbert space. In loop diagrams with two internal 
quark propagators, such as are needed for calculating mesonic bound states, one has to be 
careful about how the cut-off is imposed on both propagators if axial Ward identities are 
not to be violated. Furthermore, while a cut-off is needed to make the nonanomalous part 
of the effective action finite, cutting off the anomalous part means that, for example, the 
amplitude for vr" -^ 27 is not given correctly.^ 

An attractive alternative to the usual local NJL model is suggested by the instanton 
picture of the QCD vacuum studied by Diakonov and Petrov P,^. The interactions of 
the quarks with the instantons generate an effective four-point coupling (in the case of two 
light quark flavours). The nonlocality of this interaction provides a natural cut-off on the 
loop integrals. Moreover the nonlocality a separable form. This considerably simplifies 
the Schwinger-Dyson and Bethe-Salpeter equations in comparison with models based on 
quark-quark forces motivated by gluon exchange. [| 

In Diakonov and Petrov's original version the form factor in this interaction is given by 
the zero mode of a quark in the presence of an instanton. Here we consider more general 
choices, motivated in part by the fact that the confining nature of the QCD vacuum will 
modify the long-range behaviour of the form factor. We explicitly construct a conserved 
axial current for this type of model and demonstrate that PCAC, implicitly assumed by 
Diakonov and Petrov, does indeed hold. In addition we show that a rather natural choice of 
form factor can lead to quark confinement. We also compare our model with other nonlocal 
NJL-like models based on separable interactions. 



"^For a review of these problems, see Q. 

^See for a recent example of such a model and many further references. 



II. NONLOCAL MODEL 

The action for a generalised NJL model with a nonlocal four-point interaction and ap- 
proximate SU(2)xSU(2) chiral symmetry is 

S= d'hc %jj{x){i^ — mc)^{x) (1) 

X ^{Xi)ip{x3)^{x2)i'{x4) + ^iXi)i'y5T°'^{x3)ip{x2)i'j5^"'^{x4) 

The symmetry is explicitly broken by the current quark mass rric- For a separable interaction 
the smearing function may be written in momentum space as 

a{pi,P2,P3,PA) = |(2vr)®G'/(pi)/(p2)/(P3)/(P4) S{pi+P2-P3-Pa), (2) 

where we take the form factors f{p) to be normalised such that /(O) = 1. In the limit of a 
local interaction, f{p) = 1 for all p, this reduces to the ordinary NJL model with coupling 
strength G. 

As in the ordinary NJL model, for a sufficiently large coupling the quarks acquire a 
dynamical mass. With a nonlocal interaction this dynamical mass M{p) is momentum- 
dependent. It can be determined, in the Hartree approximation, from the gap equation: 

2 f d^k M{k)p{k) 



MW = ™. + G./(rt^yp^jj^^. (3) 

where z/ = ANcNf, Nc and Nf being the numbers of colours and flavours respectively. The 
momentum integral has been Wick rotated so that the four-momentum k is in Euclidean 
space. From now on all momenta will be Euclidean, unless explicitly stated otherwise. 
The separable nature of our interaction means that the loop integral is independent of the 
external momentum p and so this equation can be solved in exactly the same way as in the 
NJL model. 

In choosing the form factor f{p), we obviously want it to be Lorentz invariant. We 
assume that the interaction is generated by nonperturbative features of the QCD vacuum 



and so does not have the divergent short-distance behaviour of one-gluon exchange. Hence 
we want a form that falls off sufficiently fast at high momenta to keep the loop integrals finite. 
At low momenta we want it to be more regular than the instanton form used by Diakonov 
and Petrov since in coordinate space it corresponds to a function that falls off rapidly with 
distance. Such a fall-off is to be expected from the confining nature of the QCD vacuum. The 
simplest choice satisfying all of these requirements is a Gaussian function of the Euclidean 
four- momentum: 

/(p) = exp(V/A^), (4) 

where A is our cut-off parameter describing the range of the nonlocality. We use this form 
factor in the numerical calculations presented here. Other choices which also possess these 
features are likely to lead to qualitatively similar results. 

With this choice of form factor, the model does exhibit confinement of the quarks for 
large enough values of the dynamical mass. This arises because the equation p^ = M{pY = 
M(0)^exp[2p^/A^] in Minkowski space has no real solutions. Hence the quark propagator 
has no real poles and quarks do not appear as asymptotic states. The quark propagator 
does have a pair of poles with complex masses, corresponding to quarks which have a finite 
lifetime as isolated particles. This is similar to the situation in other models for confinement 
based on approximate Schwinger-Dyson equations @,0]. 

The model has three adjustable parameters: G, A and rric. We fix two of these using the 
pion mass and decay constant. The remaining free parameter can then be chosen to give 
reasonable values for the dynamical quark mass and decay constant. 

III. BETHE-SALPETER EQUATION 

The pion is constructed, as in the NJL model, by solving the corresponding Bethe- 
Salpeter equation. For a separable interaction this equation is just a geometric series. For 
example, the T-matrix in the pseudoscalar isovector channel can be written in the form 



nPl,P2,P3,P4) = f{Pl)f{P2)f{p3)f{p4)[G + GVpp(p) + G'jppip)' + ■■■] (5) 

= fip^)fip^)fip^)fip^ K_GJMpy 

where p = pi + P2 = Ps + Pa, and Jpp is given by a quark loop with two nonlocal pseudoscalar 
insertions: 



Jpp(p) =tr / 



{2ny 



r{k + \p) ., f\k-\p) 



«7 ,, , 1 y , „w, , 1 N ^7 



(6) 



1^ +\i +M{k + \p) ' 1^ -y + M{k - ip) 

The pion mass and wave function renormalisation can be determined from the position and 
residue of the pole in the propagator in this channel. For momenta in the vicinity of the 
pion pole, p^ = — m^, we write 

- = ^- +■■■. (7) 

1-GJMp) P^ + my ^> 

Assuming that the symmetry-breaking current quark mass is much smaller than the 
dynamical mass, we now expand Jpp about the chiral limit. To first order in the current 
quark mass and p^ it can be written in the simple form 

jpp(p)~l_^!!^(^^)_^ (8) 

where Mo(p) is the dynamical mass at zero current quark mass. In obtaining this we have 
made use of the gap equation and the expression for the quark condensate, 

in the chiral limit. To leading order in rric, the pion mass is given by 

ml = -J^2^Mi^). (10) 

and the wave function renormalisation by 

^1 _ z/ r d^k M{kf - M'{k)M{k)e + {M'{k)fk^ 

^ 'mW'W? (p + Moikyf ' ^^^^ 

where the prime denotes differentiation with respect to k"^. 



IV. PION DECAY CONSTANT 

In order to calculate the pion decay constant F^ for this model we need to construct an 
axial current which is conserved in the chiral limit. The usual local current, 

j^iioc) = l^7^75r>, (12) 

is not conserved in the presence of the nonlocal interaction of Eq. (1) and so it cannot be 
the symmetry current of this model. This can be seen from the results of Diakonov and 
Petrov where the pion decay constant calculated from this current fails to satisfy the 
Gell-Mann-Oakes-Renner (GOR) relation 0. In addition, the correlator for this current is 
not transverse, another signal that it is not the appropriate symmetry current. 

Various techniques are available for construction of a conserved current in the presence 
of nonlocal interactions. A rather elegant one is to make the action locally invariant by 
coupling the quarks to right- and left-handed gauge fields, A'^^ j^ [|10| . Some of these fields 
can be identified with those of the electromagnetic or weak interactions, but in general they 
need not be physical. The corresponding currents are obtained by differentiating the action 
with respect to the gauge fields and then setting those fields equal to zero. 

The locally invariant action is constructed by replacing the ordinary derivative by the 
covariant one in the kinetic term. The nonlocal interaction can be made invariant by intro- 
ducing parallel transport operators. 



WR^L{x,y) = Pexp 



1 (iz^^ 

dX—A,{z){l±^,) 
aA 



(13) 



where the integral follows a path z{X) between the appropriate right- or left-handed quark 
fields at space-time points x and y. This is analogous to the gauge-invariant point splitting 
used in deriving anomalies ||rT|. Although this method is convenient in a version of the 
model (1) with only one quark flavour, it becomes rather cumbersome for the in a model 
with two or more flavours because of the non-abelian nature of the chiral symmetry. 

Alternatively one can use a Noether-like method of construction, evaluating the diver- 
gence of the local current with the aid of the equations of motion obtained from the action 
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(1) (with rric = 0). For the axial isospin current, which we are interested in here, this gives 

^M (|^7^75T">) = - / d'^xi . . . d'^Xi a{xi,X2, X3, X4)ip{xi)Tai'y^ip{x3)ilj{x2)ip{x4) (14) 

X [6{xi — x) — 6{X2 — x) + 6{X3 — x) — 5(X4 — x)] 
+ \ I d'hci . . .d'h:ia{xi,X2,xz,X4)'il){xi)[Ta,n]'ip{x'i)'ip{x2)ni^^'il){xA) 
X [5{xi — x) — 5{x2, — x)] , 

where we have assumed that a is symmetric under permutations of its arguments. The 

r.h.s. of this equation, which arises from the variation of the nonlocal interaction under local 

axial isospin transformations, must be re-expressed as a divergence if we are to identify 

a conserved current. Ball and Ripka ||^ do this by making a derivative expansion of the 

nonlocality. 

However such an expansion is not essential: a trick suggested by the previous method is 

to write the differences of delta functions as line integrals, for example, 

/■I dz^ d 
S(xi — x) — 5(x4 — x) = / dX——— — S(z — x), (15) 

JO dX oxf^ 

where the path z{X) runs from xi to X4. Interchanging the order of differentiation with 

respect to x and the various integrals leaves the r.h.s. of (14) in the form of a total divergence 

from which the the nonlocal contribution to the symmetry current can be identified. If we 

choose, for convenience, straight line paths such as 

z{X) = {1 - X)xi + Xx4, (16) 

then the nonlocal current is 

isCni) = d^xi.-.d'^Xi / dXa{xi,X2,X3,X4) (17) 

X \il>{xi)Tai'y^ip{x3)'ijj{x2)^{x4) (^4 — Xi)6(Xxi + (1 — A)X4 — Xj 

- {^3 - x^)6(^Xx2 + (1 - A)x3) - a;) 
- ie''''''i^{xi)T^ij{x3)i^{x2)ni-f^^ix4)ix^ - x1)6[Xxi + (1 - A)a;3 -x)]. 

The full current is the sum of this and the local piece (12). 



Note that none of these procedures determines the current uniquely. The ambiguity 
shows up most clearly in the arbitrary choice of path for the line integrals but it affects 
all methods for constructing a conserved current from a nonlocal action. It arises because 
the requirement that the divergence of the current vanish only fixes the longitudinal part of 
the current; the transverse part remains undetermined. In the derivative expansion of [^] it 
corresponds to the fact that one can always add to the current combinations of field deriva- 
tives whose divergence vanishes by construction. This problem is well known in nuclear 
physics: the longitudinal components of exchange currents can be related to phenomeno- 
logical nucleon-nucleon forces, while the transverse currents require a specific model for the 
underlying meson exchanges.^ For our present purpose, calculation of F^, only the longitu- 
dinal current is needed and so our results are not affected by this ambiguity. 

The pion decay constant, —iF^p^' = (0|J5"^(0)|7r''(j9)), can be found by extracting the 
pion pole from the correlator of two axial currents. This correlator can be written 

where Jpp(j9) is given by the quark loop described in the previous section and Jpa(p) is 
given by a similar loop with one pseudoscalar and one axial current insertion. Equating the 
residues at the poles on both sides, one finds 

4F^Vp^ = z^jUpYjUp), (19) 



where p^ — --^^ 



m;. 



The contribution of the local piece of the current to Jp^ip) is straightforward and has 
the form: 



Jr 



PA(ioc)(p)-tr J ^2^^4 



d^k 



Yl^ u . 1 i . ../, . 1 n H 



1 . 5 f{k + \p)f{k-\p) 



lit +y +M{k + y) ' It -y + M{k - \p) 



(20) 



^See, for example, 



The nonlocal contribution is more involved and so we do not present the details of its 
evaluation here. Its main features can be seen from the form of the the nonlocal current 
(17) in momentum space, 



X 



f{h + P)f{h)f{h)f{h) - f{h)f{h + P)f{h)f{k 



4; 



+f{h)f{h)f{h + P)f{h) - f{h + P)f{h)f{h)f{h + P) 



+ 



where the terms involving commutators of isospin matrices have not been written out since 
they do not contribute to F.,^. Each of the terms shown in Eq. (21) involves the product 
of pseudoscalar isovector and scalar isoscalar densities. The vacuum to one-pion matrix 
element of the nonlocal current consists of similar terms, each of which is a product of two 
integrals: one similar in its Dirac structure to Jpp, the other like the self-energy in Eq. (3). 
After a certain amount of effort, the total J^x{p) can be expressed in terms of the same 
integral as in Z^^ above, and so can be written as 

J^^{p) = -2tp^Mo{0)Z-' (22) 

to lowest order in 171^. Using this in Eq. (19), we see that F^ can be related to the wave 
function renormalisation by 

F^ = Mo{OfZ;\ (23) 

The GOR relation follows at once from this and Eq. (10) for the pion mass: 

F^ml = -mS^)- (24) 

V. RESULTS AND DISCUSSION 

We fix the current quark mass nic and one combination of G and A to give a pion mass 
of 140 MeV and a pion decay constant of 93 MeV. Results for the dynamical quark mass 
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and the quark condensate are shown in the table. The variation of these quantities with the 
cut-off A is quahtatively similar to that obtained using the local NJL model with the usual 
regularisation schemes 0,^. 

Typically the dynamical quark mass in quark models is taken to be in the region 300-400 
MeV. The average current mass of the up and down quarks is believed to be about 5-10 
MeV at the momentum scale relevant to hadronic models ||13|. From the GOR relation, 



this corresponds to a quark condensate in the range (gg) = \{4'4') — —(210 MeV)^ to 
— (260 MeV)'^. From the table we see that parameter sets with cut-offs of around 1 GeV 
yield reasonable results for quark masses and the condensate. 

For large enough values of the dynamical mass, the quark propagator ceases to have 
poles at real p^ in Minkowski space, and quarks become confined. In the chiral limit, this 
happens for 

^ > 5Z. (25) 

A 2 ^ ^ 

Neglecting the effect of the small current quark mass, we see that our quarks are confined for 
those parameter sets which give Mo(0) above about 300 MeV. The critical point for confine- 
ment is thus in the middle of the region which yields reasonable values for the condensate. 
Other related models have been studied by various authors. In particular Diakonov and 
Petrov 1^ have used a nonlocal interaction with a separable form suggested by an instanton 
approach to the QCD vacuum. The expression for F^ given by Eqs. (23, 11) agrees with 
theirs 0. However those authors never explicitly calculate F-j^ from the conserved axial 
current. Instead they take the pion wave function renormalisation, and assuming that 
PCAC holds, determine Ft^ from it using Eq. (23) (the equivalent is Eq. (30) of their paper). 
When they attempt to calculate F^ using the local piece of the current alone, they find an 



expression which was originally obtained by Pagels and Stokar [Q. The incompleteness of 
this expression of Pagels and Stokar and the fact that the full symmetry current requires a 
nonlocal contribution have also been pointed out by Holdom and coworkers in the context 



of a constituent quark model with a nonlocal self-energy [|15[. In fact the model of |I5| has 
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many features in common with the one we study here, although it is not derived from an 
underlying quark-quark interaction. Ball and Ripka [Q discuss in detail the construction of 
conserved currents in such models. 

Buballa and Krewald [|I^ have also studied a model of the form (1), with a different 
choice of form factor. They too find quark confinement but their form factor is chosen to 
avoid the appearance of complex poles in the quark propagator. The propagator in their 
model does however contain some unusual analytic properties, in the form of additional cuts 
arising from the square roots in the form factor. The appearance of extra singularities in the 
quark propagator seems unavoidable in models for confinement based on Schwinger-Dyson 
equations @,0]. Moreover the model of Ref. |T^ still requires an explicit cut-off since the 



form factor fails to regulate the momentum integrals in loops. 



Another related model has been studied by the Rostock group [0. It is based on 
an instantaneous interaction with a separable dependence on the three-momenta. That 
model is of course not covariant and so, for example, there will be different pion decay 
constants associated with the space and time components of the axial current. Schmidt et 
al. have calculated F^ from the time component; since their interaction is instantaneous, 
that component of the current is purely local and the extra terms described here do not 
appear. 
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TABLE 



Mo(0) 



A 



-{qq) 



rrir 



200 
250 
300 
350 
400 
450 
500 



2222 
1414 

1075 
887 
768 
684 
623 



330 
257 
222 
200 
185 
173 
165 



2.4 
5.0 
7.7 
10.5 
13.3 
16.2 
18.9 



TABLE I. Results for dynamical quark mass and quark condensate, with Ft^ = 93 MeV and 
?TT-,r = 140 MeV. All quantities are in MeV. 
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